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ABSTRACT 

We investigate the evolution of circumstellar discs around T Tauri stars that are 
ejected from small stellar clusters within molecular clouds. In particular, we study 
how the interaction that leads to ejection may hasten the transition between Classical 
and Weak-lined T Tauri status. In our models, ejections of T Tauri stars at velocities 
of 3-10 km/s truncate the accretion disc at radii between 1 and 10 a.u., reducing the 
viscous evolution time of the disc so that accretion rapidly ceases. The observational 
appearance of the resulting systems is then dependent on the presence or absence of 
a stellar magnetic field. For non-magnetic stars we find that a near-infra red excess 
should persist due to reprocessing of stellar radiation, but that this is greatly dimin- 
ished for magnetic T Tauri stars by the presence of a magnetosphere extending to 
corotation. In either case, there is a period when ejected stars should appear as non- 
accreting systems with detectable circumstellar material at wavelengths of 5 microns 
and beyond. We discuss the implications of these results for models in which ejected 
stars contribute to the halo of pre-main-sequence objects discovered from ROSAT 
observations of star forming regions and the All-Sky Survey. 
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1 INTRODUCTION 

The ROSAT All-Sky Survey (RASS) has identified large 
numbers of candidate Weak-lined T Tauri stars (WTTS) in 
and around nearby star- forming regions, including Taurus- 
Auriga (Neuhauser et al. 1995a), Chamaeleon (Alcala et al. 
1995) and Orion (Alcala 1994; Sterzik et al. 1995). Spec- 
troscopic follow-up of subsamples of these candidates finds 
strong Li absorption in around 60% of surveyed sources (Al- 
cala et al. 1995, 1996), and in some cases provides kinematic 
evidence that the stars newly identified from the RASS 
have radial velocities consistent with the previously-known 
WTTS population (Neuhauser et al. 1995b). In most case, 
however, such evidence is lacking, and indeed the most strik- 
ing observation is the large spatial extent of the candidate 
WTTS, many of which lie far (10° or more) from known star- 
forming regions (Montmerle & Casanova 1995; Neuhauser et 
al. 1995b; Sterzik et al. 1995). 

The distances of the newly identified sources are mostly 
highly uncertain, making age determinations from the po- 
sition of stars in the Hertzprung-Russell diagram suspect. 
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Montmerle & Casanova (1995) quote very young ages (~ 10 5 
yr) for some sources relatively close (a few degrees) to 
Chamaeleon, while Neuhauser et al. (1995b) give upper lim- 
its of 25 Myr for a sample of WTTS south of Taurus based 
on Li abundances. For a larger sample of ROSAT discovered 
sources, Magazzu et al. (1996) find an overabundance of Li 
as compared to the Pleiades for a range of spectral types, 
including those (eg. K and later) for which the Li width be- 
comes important as an age indicator. These results suggest 
that the selection criteria for the larger samples are finding a 
new population of WTTS younger than the Pleiades, though 
contamination from much older (10 8 yr) stars may also be 
important (Hartmann, private communication), especially 
for the purpose of statistical studies. 

Here, we note that if any significant fraction of the 
newly identified WTTS are young (i.e. have ages in the 1- 
10 Myr range common for previously known WTTS), then 
their origin so far from known sites of star formation poses 
two severe problems. Firstly, how did these sources arrive 
at their current location, and secondly, why is a correspond- 
ing population of dispersed Classical T Tauri stars not also 
observed. One possibility is that the number of stars form- 
ing in small isolated groups has previously been underes- 
timated, as suggested by Feigelson (1996) on the basis of 
simple kinematic models for the distribution of the dispersed 
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RASS stars. Alternatively, the halo sources might have been 
formed within known regions of star formation activity, be- 
fore being ejected by dynamical interactions. In this latter 
scenario the number and observed distribution of the out- 
liers implies that a reasonable fraction of all low-mass stars 
must be ejected from molecular clouds at relatively high ve- 
locities of several kms -1 or greater. 

The decay of small stellar clusters has been considered 
as a contributor to the halo of ROSAT sources by Sterzik 
& Durisen (1995), and in the context of binary formation 
by McDonald & Clarke (1995). In these calculations decay 
of the cluster and ejection of most of its members occurs 
rapidly - within a few tens of crossing times (a few x 10 5 
- 10 yr for the McDonald & Clarke simulations, and ear- 
lier for the denser clusters modelled by Sterzik & Durisen). 
Although the lifetime of circumstellar discs varies widely 
(Strom 1995), at such early epochs most of the ejected stars 
would be expected to be Classical T Tauri stars (CTTS), 
surrounded by circumstellar discs. Since CTTS are generally 
well-localised in dark clouds (eg. Alcala et al. 1995), a viable 
ejection model for the halo ROSAT sources requires that the 
ejected CTTS are converted into exclusively Weak-lined sys- 
tems (where there is little or no evidence for circumstellar 
material) on a timescale that is short compared to the ages 
of the halo sources. 

In this paper, we discuss the effects of an encounter lead- 
ing to ejection on the accretion disc of a CTTS. The ejection 
of stars with velocities of a few kms -1 implies interactions 
that are close enough to severely curtail the extent of the 
accretion disc, but not so close as to destroy it entirely. The 
weakened disc will then evolve more rapidly because of the 
reduced viscous timescale at its outer edge, so that the rate 
of accretion onto the star falls rapidly. However significant K 
excesses could still arise from reprocessing of the stellar lu- 
minosity by a passive disc, and we show that the strength of 
this contribution to the spectral energy distribution depends 
critically on the magnetic field of the star. For weakly mag- 
netic stars (where the disc extends to the stellar equator) 
the passive flux from the disc will persist long after active 
accretion has ceased, whereas a strong stellar field can dis- 
rupt the hot inner disc regions and reduce the reprocessing 
flux coincident with the decline in the mass accretion rate. 
Magnetic fields of the required strength have been directly 
detected in several WTTS (Basri, Marcy & Valenta 1992; 
Guenther & Emerson 1995, 1996), and many indirect lines 
of evidence suggest that they are common in both WTTS 
and CTTS (see, eg, the review by Hartmann 1994). 

The plan of this paper is as follows. In Section 2, we 
estimate the accretion rates below which magnetic and non- 
magnetic systems would appear as WTTS on the basis of 
their near-infrared colours and Ha equivalent width. In Sec- 
tion 3 we use these estimates to calculate the timescale for 
the CTTS — ► WTTS transition, and compute the stellar den- 
sities required to eject a large number of T Tauri stars from 
their parent clouds. In Section 4 we verify these timescale 
estimates by directly computing the time- dependent evolu- 
tion of the star-disc system for a typical choice of param- 
eters. Section 5 summarizes our findings and discusses the 
observational consequences of the model. 




log (accretion rate / M yr ') 

Figure 1. K excesses over the stellar photosphere for a variety of 
accretion disc models described in the text, as a function of the 
steady-state accretion rate through the disc, (a) A non-magnetic 
disc extending to the stellar surface, (b) A magnetically disrupted 
disc around a star of period P* = 4 dy with P m = R c (lower 
curve) and P m = 0.7 X R c (upper curve), (c) A magnetically 
disrupted disc around a star of period P* = 8 dy with P m = R c 
(lower curve) and P m = 0.7 X R c (upper curve). The curves for 
magnetic discs with P m = 0.7 X R c do not extend to very low M, 
since at low accretion rates the magnetosphere must lie at larger 
radii. 

2 MODELS FOR T TAURI ACCRETION DISCS 

The parameters and mechanisms that control the transi- 
tion of a star from Classical to Weak-Line status are largely 
unknown, but it is natural to assume that the mass accre- 
tion rate through the disc M plays an important role. With 
this assumption, we define three mass fluxes; the accretion 
rate for a typical CTTS, Mqtts, the accretion rate below 
which a star appears as a WTTS in a non-magnetic model 
Mweak.noB, and the corresponding quantity in a magneto- 
spheric model Mwoak.B- A variety of observational charac- 
teristics distinguish WTTS from CTTS (see, eg. the review 
by Bertout 1989), but low equivalent width of Ha ( < 10A) 
and lack of an infra-red excess at K (2.2 fim) are common 
identifying features. We adopt these conditions for deciding 
when model systems have become WTTS, and estimate the 
values of Af wca k, n oB and M woa k,B- These quantities are then 
used to estimate the timescale for the CTTS to WTTS con- 
version in magnetic and non-magnetic systems that suffer 
an ejection interaction. 

Figure 1 shows the excess flux AK at 2.2 pirn for a range 
of discs around a star of mass 1 Mq and radius 2 Rq. The 
models include the heating of the disc from accretion, re- 
processing of stellar radiation, and where appropriate work 
done by magnetic torques (assuming a stellar dipole field 
of 1 kG). The models are calculated using the formalism 
outlined in the Appendix. No contribution from a boundary 
layer or magnetically funnelled accretion shock is included as 
these components of the total system luminosity are not ex- 
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pected to be strong at K (note that in magnetospheric mod- 
els the infalling material is unlikely to be optically thick at 
K even for M = Mctts , and will certainly have tk <S 1 for 
the lower accretion rates of interest here). For the magnetic 
models we have assumed that the disc is disrupted inside 
some magnetospheric radius P m , which is taken either to 
be at corotation R c (where the Keplerian disc has angular 
velocity Q = SI*), or at 0.7 x R c . P m « R c is expected on 
theoretical grounds at low M (eg Armitage & Clarke 1996; 
see also Wang 1996), while fitting of magnetic disc mod- 
els to the spectral energy distributions of CTTS typically 
gives P m « 0.7 x R c (Kenyon, Yi & Hartmann 1996). Ad- 
ditional angular momentum loss via stellar winds (e.g. Tout 
& Pringle 1992) would be expected to reduce P m below the 
values seen in models that assume magnetic disc linkage as 
the only braking mechanism. The stellar rotation rate P, 
is taken to be 4 dy or 8 dy, these values being typical for 
WTTS and CTTS respectively (Bouvier et al. 1995). 

For the non-magnetic disc model, reprocessing of the 
stellar luminosity ensures that an easily detectable K excess 
is found for all accretion rates, provided only that the disc re- 
mains optically thick at K. This will be the case down to very 
low accretion rates. To order of magnitude, an optical depth 
tk = 1 requires a disc column density E = 0.1 g/cm , which 
is at least 4 orders of magnitude below the surface density 
seen in CTTS disc models accreting at 10~ 7 Mq yr -1 . The 
accretion rate corresponding to S = 0.1 depends on the disc 
viscosity v via M ~ 37r^E, where v is typically an increas- 
ing function of E. We conclude that M wea k,noB is likely to 
be more than 4 orders of magnitude below Mctts , and that 
a disc around a non-magnetic T Tauri star would remain 
visible via its K excess long after significant accretion (evi- 
denced by high equivalent width of Ha) had ceased. 

For the magnetic models, the behaviour of AK as a 
function of M is very different. If P m = 0.7 x R c , then signif- 
icant K excesses of ~ 1 mag are obtained at accretion rates 
of 10 -7 Mgyr -1 and above for P* = 4 dy, with smaller but 
still readily detectable AK values for the longer rotation pe- 
riod. However, once the accretion rate has fallen sufficiently 
low that the magnetosphere lies at corotation, the models 
predict that AK should be low. The dependence of P m on 
M is (Clarke et al. 1995), 

Rm = (^pLf\ n> 

\My/GMlJ 

Thus, taking as observationally determined the result that 
Rm « 0.7 x R c in the Classical T Tauri phase, we find that 
the magnetosphere should lie at corotation once M < 0.3 x 
Mctts- From the Figure, we then conclude that T Tauri 
systems rotating at P* = 8 dy would appear as WTTS on 
the basis of their K excess below an accretion rate of a few 
x 10~ 8 M Q yr _1 . For the shorter rotation period of 4 dy, AK 
falls below 0.1 magnitudes at approximately the same M. 

These results imply that judged by the near infra-red 
excess, a magnetic CTTS would appear indistinguishable 
from a WTTS at an accretion rate of ~ 10~ 8 Mq yr -1 
or below. The equivalent width of Ha in WTTS shows 
considerable scatter (Alcala et al. 1993), so accretion at 
approximately this level probably cannot be excluded on 
the basis of the Ha measurements either. Hence, we esti- 
mate that M woa k,B ~ 10~ 8 M© yr -1 , and conclude that 



Mveak.B S> M woa k, n oB • This critical accretion rate is then 
used in the following Section to estimate the timescale on 
which interactions leading to ejection from a cluster will 
reduce M below M wcakjB , and hence convert a CTTS to 
Weak-lined status. 



3 TIMESCALE FOR THE CTTS -> WTTS 
TRANSITION 

The decay of small stellar clusters has been studied by a 
number of authors using N-body simulations (van Albada 
1968; Sterzik & Durisen 1995), and techniques that attempt 
to incorporate additionally the effects of discs (McDonald & 
Clarke 1995) or gas (Bonnell et al. 1996). Such clusters are 
found to dissolve within a few tens of crossing times, leading 
to the ejection of most of the members and the formation 
of one central binary. The median velocities of the escapers 
Reject follow the scalings established for the 3-body problem, 

fe ^°' 5X ((S)) 1/2 ' (2) 

where m c j oct are the masses of the ejected stars, and |Po| oc 
M?/P c is the total energy of the system with mass M c 
and scale length P c (Valtonen & Mikkola 1992; Sterzik & 
Durisen 1995). For a small cluster containing N stars of mass 
M„, the above expression implies, 

This broadly agrees with the numerical results suggesting 
that stellar separations in the 300-1000 a.u. range are re- 
quired to generate significant numbers of escapers with inject 
greater than a few kms" 1 . The required clusters are thus 
vastly denser than the stellar density in regions such as 
Taurus- Auriga, even when account is given to the observed 
degree of subclustering (Gomez et al. 1993). However such 
clusters are also short lived, typically < 10 5 yr (Sterzik & 
Durisen 1995), and hence the observed state of Taurus can- 
not necessarily be taken as evidence that stars were not born 
in much more compact groups. Indeed the dissolution of 
such initial clusters might plausibly lead within a few Myr 
to both high velocity escapers and dilute aggregates (com- 
posed of stars ejected at only ~ a kms -1 ) similar to those 
observed by Gomez et al. (1993). Detailed population syn- 
thesis would be required to address this possibility in a more 
quantitative fashion. As discussed by McDonald & Clarke 
(1995) and Sterzik & Durisen (1995), the formation of stars 
within such transitory clusters is consistent both with some 
theoretical models of star formation (Boss 1993), and the 
observed high fraction of binaries in the pre-main-sequence 
(eg. Ghez 1996). 

For a dynamical ejection model to be relevant to the 
formation of the RASS halo of WTTS, u e ject must be a few 
kms -1 or greater. Here we assume that conditions in the 
initial stellar clusters are conducive to producing such ve- 
locities (i.e. similar to those used as initial conditions by 
Sterzik & Durisen 1995), and calculate the effects of those 
ejections on circumstellar discs. 

To relate v c j cct to the distance of closest approach in 
the encounter that led to the ejection, P per i, we make use 
of prior work on the scattering of single stars by binaries 
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(for details see, e.g. Davies 1995). The energy lost from the 
binary AE is related to the binding energy of the binary, 
Ebin, via 



AE < 0.4 x _Bbi„ = 0.4 x 



GM 2 

2dbin 



(4) 



where dbin is the separation of the binary before the en- 
counter. Making use of the further result that dbin ~ -Rperi, 
we find 



Wcjcct ~ 15 km/s x ( Y^~) 



-1/2 



(5) 



for stars of solar mass. Ejection velocities of w c j C ct = 3 — 
10 km/s, then imply i? pcr i ~ 2 — 25 a.u.. Star-disk encoun- 
ters have been extensively modelled (Hall, Clarke & Pringle 
1996; Clarke & Pringle 1993), and these simulations suggest 
that the accretion disc following the encounter has an outer 
edge at ~ R pcl i/2. The ejected stars would thus be expected 
to harbour remnant discs with outer radii R ou t ~ 10 a.u. im- 
mediately following the interaction. 

The weakened disc will then evolve faster as a conse- 
quence of the reduced viscous timescale t u at R out . The vis- 
cous timescale is given by (e.g. Pringle 1981), 



*"~ an (h) 



(6) 



where a is the Shakura-Sunyaev viscosity parameter, H is 
the disc scale height, and all quantities are evaluated at ra- 
dius R. For stars of solar mass, 

where a is the Shakura-Sunyaev viscosity parameter appro- 
priate to the outer regions of the disc following the interac- 
tion. 

Thus provided that ctR out Z 10 -3 (as is likely since even 
a = 10 -3 implies worryingly large disc masses at radii of 
several a.u.) and (H/R) > 0.1, the viscous timescale of the 
truncated disc will be short compared to the typical disc 
lifetime in CTTS. For example at 5 a.u., which corresponds 
to the outer edge of the disc in a star ejected at 5 kms -1 , 
t v ~ 2 x 10 5 yr with the fiducial parameters given above. The 
same conclusion follows by noting that unless a or (H/R) 
are strongly increasing functions of R in the outer regions 
of the disc, the radial dependence of t v is approximately as 
fi _1 oc R A/2 . A typical T Tauri disc has a lifetime of a few 
Myr, which implies t v ~ 1 Myr. Reducing R ou t by a modest 
factor of ~ 5 then ensures that immediately following an 
encounter (before the disc has had time to expand viscously) 
the viscous timescale will be reduced by around an order of 
magnitude. Estimates of T Tauri disc sizes are generally 
~ 10 2 a.u. or greater (McCaughrean & O'Dell 1996), and 
hence encounters at radii of ~ 10 a.u. are likely to meet this 
condition. 

The viscous timescale can be regarded as a characteris- 
tic time for the decay of the disc surface density and accre- 
tion rate. As Mctts ~ 10 -7 Mq yr -1 , only a few viscous 
times will be required before M through the truncated disc 
falls below M wca k,B- Thus within a few times t v (equation 
Q), the initially Classical T Tauri star will lose the AK and 
Ha signatures of a CTTS, and will appear as a Weak-Lined 
system. In this model circumstellar material will still exist 



around the star for a further, perhaps lengthy, period, but 
the low accretion rate and the hole created by the stellar 
magnetosphere renders it undetectable in the near infra-red. 

We emphasise that the same conclusion does not follow 
for non-magnetic disc models. In this scenario, the interac- 
tion reduces t v exactly as for the magnetic case, but a large 
number of viscous times are still necessary before M falls be- 
low M W eak,noB. Before the accretion rate falls this low, the 
disc will have had time to re-expand to its original extent, 
restoring t v to a much larger value. As a consequence, the 
reduction in the duration of the CTTS phase will be much 
less marked, and we would expect the ejected systems to dis- 
play infra-red excesses from passive discs without significant 
active accretion. 



4 NUMERICAL SIMULATIONS 

To verify the estimates derived in Section 2, we have calcu- 
lated explicitly the evolution of the coupled star-disc system 
following an encounter that destroys the outer accretion disc. 
The response of the disc to magnetic and viscous torques 
is described using the evolution equation given by Livio & 
Pringle (1992), 

aS _ 3__d_ 

dt RdR [" dR x 

1 d /fi-fi„ B 2 Z R 5 ' 2 



R^JL^R 1 ' 2 ) 



+ - 



RdR \ n -u^GMl 



(8) 



which we solve together with equations for the spin period 
and magnetic field of the star. The stellar model used is a 
0.9985 Mq model calculated using the most recent version 
of the Eggleton code (Pols et al. 1995), as modified for pre- 
main-sequence evolution by Tout (private communication). 
The model and numerical details are described elsewhere 
(Armitage & Clarke 1996). 

Figure 2 shows the time-dependence of the disc mass, 
accretion rate, stellar spin period and K excess for two 
models; one in which the disc beyond 5 a.u. is lost at 
t = 2 x 10 yr (representing the effect of an ejection in- 
teraction), and a 'control' model in which the disc evolves 
undisturbed. We also show the evolution of AK for a non- 
magnetic model subject to the same disruptive encounter. 
The main parameters of the simulations are; an outer bound- 
ary condition of E = at 50 a.u., a stellar magnetic field 
B, = 1500 G x (P„/4dy), an initial disc mass of 0.3 Mq, 
and an initial accretion rate of 1.5 x 10~ 7 M©yr _1 . The vis- 
cosity prescription assumes an a in the inner disc of 10 -3 , 
and gives an evolutionary timescale for the unperturbed disc 
of ~ 10 6 yr, consistent with the inferred lifetime of CTTS 
discs (Simon & Prato 1995). The K excess AK is calculated 
using identical models to those of the previous Section, ex- 
cept that the assumption that the disc is in a steady state 
is relaxed and we use instead the computed surface density 
to obtain the active accretion component. 

Following the encounter, the accretion rate drops 
rapidly, reflecting the reduced viscous timescale at the new 
outer edge (t u decreases by a factor ~ 10 following the in- 
teraction). At late times, when the weakened disc has re- 
expanded to its previous size, the accretion rate is a factor 
of ~ 6 times lower than in the model that evolves undis- 
turbed, or approximately (1—2) x 10~ 9 M Q yr _1 at t = 2 Myr 
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Figure 2. Comparison of the evolution of an undisturbed star-disc system (dashed line) with one experiencing an encounter at t = 
2 X 10 5 yr that removes all mass outward of 5 a.u. (solid curve). Both models assume a disc disrupted at small radii by a stellar 
magnetosphere. The panels show the time-dependence of the monochromatic flux excess AK at 2.2 (im, accretion rate, and stellar spin 
period. The K excess for a non-magnetic system undergoing the same encounter as the magnetic model is also shown (dotted line). 



for these models. The rotation period of the 'ejected' star 
is marginally increased as compared to the control model, 
because the spin-up torque from accretion is reduced while 
the spin-down torque from magnetic braking is unaffected. 
The K excess drops rapidly due to a combination of lower 
M, increased ratio of R m /R c (which is ~ 1 in the weak- 
ened disc simulation), and increased R c due to the slower 
rotation. AK falls below 0.1 magnitudes within 0.5 Myr, 
and is negligible at the end of the simulation. This reduc- 
tion in AK as compared to the control run (which is itself 
still magnetically disrupted at small R) amounts to a factor 
of ~ 10, in agreement with the estimates presented previ- 
ously. The run in which a non-magnetic disc is subjected to 
an encounter also follows the expected behaviour. Immedi- 
ately after the interaction, AK drops as for the magnetic 
case due to the reduced M. However, at late times the K 
excess remains strong due to the reprocessing luminosity of 
the undisrupted disc at small radii. 

The extent to which the length of the CTTS phase is 
reduced by encounters is shown as Figure 3. A series of cal- 
culations with the parameters given above were run, varying 
only the outer radius at which the disc was truncated by the 
interaction. The time subsequent to the interaction required 
for M to fall to 10 -8 MQyr -1 is plotted as a function of 



Wcjcct, which is related to R peT i via equation (|5[). Also plot- 
ted is the time required for AK to fall below 0.1 magnitudes. 

From the Figure, it can be seen that encounters that 
produce high velocity escapers lead to a rapid decline in 
M and AK on a timescale of at most a few x 10 5 yr. The 
most destructive encounters considered here (at 10 kms - ) 
reduce the accretion rate by an order of magnitude in 10 yr. 
At lower velocities (ti e ject ~ 2 kms -1 ), the influence of the 
interactions tapers off as the mass and viscous time of the 
remnant disc approach the values of the undisturbed system. 
For the viscosity and initial conditions (primarily, the choice 
of initial disc mass) used here, this occurs at « e ject 5; 2 kms -1 . 

The form of the curves shown in Figure 3 reflects the re- 
lationship between the viscous timescae of the undisrupted 
disc and that of the truncated disc following the encounter 
(equation [?]). There are little in the way of observational 
constraints on the disc viscosity at large radius (essentially, 
just the disc lifetime which is expected to be 0(t v ) for the 
disc), and thus there is bound to be considerable uncertainty 
inherent in these timescale estimates. Nonetheless, the re- 
sults presented here indicate that encounters close enough 
to eject stars at velocities, u c j C ct ~ 5 kms -1 , should be dis- 
ruptive enough of the disc to rapidly cut-off accretion, and, if 
the star possesses a magnetosphere extending to a few stel- 
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Figure 3. Effect of the encounter on the accretion rate and K excess, as a function of the ejection velocity t> e j oc t (or, cquivalently, the 
pcriastron separation), (a) Time after encounter required for the accretion rate to fall below 10 — 8 MQ-yr -1 . (b) Time required for the K 
excess to drop below 0.1 magnitudes. In both panels the period required for the undisruptcd disc model to reach these limits is shown 
as a dashed line. 



lar radii, additionally lead to a sharp drop in the near-IR 
excess. These high velocities are required if any significant 
fraction of the most dispersed RASS sources are relatively 
young (< 10 Myr) and originated as escapers from known 
molecular clouds (Sterzik et al 1995). Less destructive en- 
counters, with periastron separations of 20-30 a.u. and ejec- 
tion velocities of ~ 3 kms -1 , are predicted to reduce the 
CTTS phase by modest factors of around 2. 



5 SUMMARY 

In this paper, we have investigated the evolution of discs 
around Classical T Tauri stars which are subjected to dis- 
ruptive encounters within the environment of small clusters. 
We find that interactions that are able to eject stars from the 
cluster at velocities of a few kins -1 or greater lead to trun- 
cation of the accretion disc, and greatly accelerate its evo- 
lution. The rate of accretion through the disc drops rapidly 
for several viscous times of the remnant disc, and falls to 
very low levels within, typically, a few x 10 5 yr. Since ejec- 
tions at these velocities from small clusters happen at an 
early epoch, the curtailment of further infall onto the outer 
regions of discs in escaping stars would also tend to reduce 
the accretion rate at later times. 

The further evolution of the discs in ejected systems 
then depends on the strength of the stellar magnetic field. 
For systems with weak magnetic fields, we predict a period 
in which a passive reprocessing disc generates prominent 
infra-red excesses at all wavelengths. Alternatively, for stars 
with a strong, ordered, magnetic field, the stellar magneto- 
sphere should rapidly disrupt the inner regions of the disc 
out to close to the corotation radius. As magnetic fields of 



the required strength are inferred from other observations, 
we regard this magnetic scenario as the more probable. Ob- 
servationally such systems would show negligible infra-red 
excesses at K, but would still possess infra-red emission at 
wavelengths A > 5 fim. Such transition systems would have 
stellar rotation periods comparable to the Classical T Tauri 
stars, and detectable emission at mm wavelengths that is, 
however, depleted by the reduction of the disc mass as a 
result of the encounter. 

In the context of models wherein the Li-rich halo RASS 
sources originate as escapers from young clusters (Sterzik & 
Durisen 1995), the present study implies that the lack of a 
dispersed CTTS population is not a problem, even if some 
fraction of the dispersed WTTS turn out to be young (~ 1 
Myr). The interactions leading to ejection at such high veloc- 
ities necessarily cut-off accretion within less than this time. 
However there should be systems (preferentially the younger 
ones), where residual circumstellar material exists, and this 
may be detectable at longer infra-red and mm wavelengths. 
Measurements of the proper motions and determinations of 
the ages of these RASS sources should provide a decisive 
test of such scenarios. 

This model, in common with previous work (Clarke et 
al. 1995), predicts stars with properties intermediate be- 
tween Classical and Weak-Lined T Tauri systems. In par- 
ticular, if the disc is cleared at the end of the CTTS phase 
'inside-out', for example by winds or a stellar magneto- 
sphere, it is hard to avoid a period in which the system 
is 'weak' on small scales but 'strong' on large scales. Such 
systems should have detectable mid infra-red and mm fluxes 
without corresponding K or Ha emission. If such systems 
are not observed, that would be strong evidence in favour 
either of an unexpectedly short disc clearing timescale, or 
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of disc dissipation via processes unconnected to the central 
star. 
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APPENDIX: CALCULATING DISC SPECTRAL 
ENERGY DISTRIBUTIONS 

The spectral energy distributions for the magnetic disc mod- 
els are calculated assuming that annuli in the disc radiate 
as blackbodies with effective temperature T c , where 



t* + t* + t£, 



0) 



and the terms on the right-hand side represent respectively 
the contributions from active accretion, reprocessing of the 
stellar luminosity, and work done on the disc by the magnetic 
field. For a disc in which the energy from viscous dissipation 
is radiated locally 



4 



(10) 



where a is the Stefan-Boltzmann constant. Evaluation of the 
right-hand side of this expression does not require knowledge 
of the specific form of u, but for a magnetic disc model it 
does depend on the assumed magnetic field configuration. 
In Section 2 we compute vT, from the steady-state solution 
of the disc evolution equation with boundary conditions 
i^E = at R = R m and vT, = M /(3n) at large radius (the 
standard non- magnetic disc expression), yielding 
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(11) 



with = (B,Rl) 2 /(Vity/GM,). The main assumption re- 
quired to derive this expression is in the form of the radial 
dependence of the magnetic torque, and is described in de- 
tail in a previous paper (Armitage & Clarke 1996). Note that 
this solution is not in general physically reasonable (yH > 0) 
for all choices of i? m , reflecting the fact that for large B* 
and/or small M steady solutions with i? m <C R c do not 
exist. In Section 3 we use vY* calculated directly from the 
time-dependent solution of equation (j§|). 

For the component due to reprocessing of stellar radi- 
ation, we utilise the results of Kenyon & Hartmann (1987) 
and Adams & Shu (1986). These authors obtain 



(1- 



(12) 



for the case assumed here of a flat reprocessing disc. 

The magnetic torques arising from field lines linking the 
star and the disc will also do work on the disc material, and 
some fraction of this may be thermalized and heat the disc 
material. Assuming all the available energy goes into heating 
the disc, 
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The actual fraction of the available energy that goes into 
reheating the disc material is uncertain, and so for the pur- 
poses of this work we adopt the conservative option of choos- 
ing the maximum value. Even then, this term is typically 
smaller than T p , and so the uncertainty in its value does not 
seriously change the estimate of the disc temperature at low 
accretion rates. 
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